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The two-dimensional (2D) Ising model (2DIM) is a one of only few nontrivial solvable 
models in statistical mechanics [1-18]. In its original formulation, this is a discrete-spin 
lattice model for order-disorder transition [1-7]. A remarkable feature, however, is that 
2DIM admits as well a field-theoretical interpretation in terms of free Majorana fermions 
on a lattice [8-13,16,17]. This significantly simplifies the analytics of the 2DIM and pro- 
vides a new insight into the physical nature of the problem. In this contribution, we 
shortly comment on some aspects of a simple non-combinatorial fermionic approach to 
2DIM based on the use of the integrals over the anticommuting (Grassmann) variables 
and the mirror-ordered factorization ideas for the 2DIM density matrix. A more detailed 
discussion on related subjects is also given in [17]. Let us start with a generalized for- 
mulation of the 2D Ising model, assuming arbitrary inhomogeneous distribution of the 
bond coupling parameters over a rectangular lattice net. The Ising spins, a mn = ±1, 
are disposed at the lattice sites, mn, with m,n = 1, ...,L running in the horizontal and 
vertical directions, respectively. At final stages, N = L 2 oo. The hamiltonian is: 
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where 6^ = (3J£°$ are dimensionless bond coupling parameters, (3 = 1/kT is the inverse 
temperature. The free boundary is assumed: <JL+in = 0, u m L+i = 0. The partition 
function is: Z = E exp (—/3H), where the sum is taken over the all possible spin config- 
urations provided by a mn = ±1 at each site. Noting an identity for a typical Boltzmann 
weight, exp (baa') = cosh b + aa' sinhfr, which readily follows from [aa') 2 = 1, we come 
to the reduced partition function: 
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where t^) = tanh b^^ , and Sp^ stands for a properly normalized spin averaging, such 
that Sp(l) = 1, Sp (a mn ) = at each site. We intend to convert Q into a Gaussian 
fermionic integral. The starting point is the fermionic factorization of the local Boltz- 
mann weights from (2) [11,17]. For the whole lattice, we introduce a set of the totally 
anticommuting Grassmann variables, a mn , a^ n , b mn , b^ n , and write: 
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In a conventional notation, the bond Boltzmann weights are now presented as A mn A* l+ln 
and -B mn -B^ n +i, where the separable factors (to be called shortly Grassmann factors) are 
to be identified from (3). The identities (3) readily follows from the elementary rules of 
fermionic integration. 2 At the next stage, we keep in mind to group together, over the 
whole lattice, the four factors with the same spin (the same index mn) and to sum over 
a mn = ±1 in each group of factors, independently, thus passing to a purely fermionic 
expression for Q. The four Grassmann factors with the same spin are: 
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These factors are coming by factorization of the four different bonds adjacent to a given 
mn site. The problem is, however, that the separable Grassmann factors like (4) are, in 
general, neither commuting nor anticommuting with each other. So, a special ordering 
of factors in their global products is needed in order the elimination of spin variables 
be really possible [11, 17]. 3 In two dimensions this problem is solvable by means of the 
mirror-ordered arrangements of Grassmann factors in the total products of weights (3) 
forming the density matrix in (2) [11,17]. Thus, we write: 
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where Sp( a6 ) stands for the diagonal Gaussian averaging arising in (3). At any given 
fixed n, we find here the two ordered m-products of Grassmann factors, which are then 
multiplied over n (the modes of ordering are shown by arrows). The averaging over 
a mn = ±1 is to be performed at the junction of the two m-ordered products in (5), 
for given n. At each step, we average over a rnn = ±1 the product A^ n S^ n A mn S mn 
at the junction. This results an even fermionic polynomial, equivalent to the Gaussian 
exponential factor just corresponding to the local term (the mn-term) in the second line 
of (6). The equivalence of the polynomial to the exponential can be checked either by a 
direct calculation, taking into account the nilpotent properties of fermions, or making use 
of the identities like (11), see below. By a repeating use of the same procedure, taking at 
each stage the commuting exponential factor away from the junction, we finally eliminate 
all spin variables in (5) [11,17]. The partition then appears in the form: 



Q = J II II da mn da mn db mn db mn eX P { J2 J2 [ a mn^mn + KnKn + 
m=l n=l m=l n=l 

_i_/(l)/(2) n * U* i fi(l)„« ++(2)/,* W„ A_h ) + n h U (6) 
' v mn v mn "m-ln^nm-l ' \ v mn u 'm—ln ' v mn u mn~\) V^mn ' u mn) ' u 'mn u mn J J i V u / 



Let us remember that Grassmann variables are the purely anticommuting fermionic symbols. Given 
a set of Grassmann variables, a\, ai, . . . , ajv, we have: + cijdi — , a"j = 0. The rules of integration 
for one variable are: J daj ■ aj = 1 , J daj ■ 1 = [8]. The Gaussian fermionic integrals of any kind are 
expressible in terms of the determinant (Pfafnan) of the correspondent matrix [8] . 

3 For a more detailed comment on this subject see [11,17]. The construction of the mirror-ordered 
representation (5) for the density matrix is the clue point of the solution of the 2D Ising model in a zero 
magnetic field within given approach [11,12,16,17]. The ordering problem is an obstacle, however, to 
solve the 3D Ising model. A comment on the 2D Ising model in a nonzero magnetic field is given below. 
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which is a Gaussian fermionic integral. Equivalently, the 2DIM is reformulated as a theory 
of free fermions on a lattice. The above expression for Q is exact, assuming free boundary 
conditions for fermions: ag„ = , b^ = 0. For further transformations of the integral (6) 
(reduction to two variables per site; continuum limit; effects of disorder) see [16,17]. 

Let us now consider the homogeneous rectangular lattice, t£l = t 1: t^l = t 2 . The 
integral (6) can be explicitly evaluated in this case by passing to the momentum space 
for fermions. In the momentum space, the integral becomes: 
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where a„ n , a* , b„ n , b* are the new variables of the integration introduced by the standard 
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Fourier substitution with periodic boundary conditions (change of boundary conditions 
can be viewed as a boundary approximation inessential in the limit of infinite lattice). 
After a proper symmetrization of the fermionic sum (action) from (7) with respect to the 
p, q <-> L—p, L—q conjugation, the integral decouples into a product of the elementary low- 
dimensional integral factors, Q 2 . These integral factors, Qp q , can be readily evaluated by 
the standard rules of fermionic integration [12,17]. Finally, we obtain an explicit solution 
for the squared partition function in the form: 

Q 2 = nn [(l+t 1 2 )(l + t 2 2 )-2t 1 (l-t 2 2 )cos^-2t 2 (l-t 1 2 )cos^l. (8) 
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This is the exact solution for Q 2 of the 2D Ising model in the limit L 2 — > oo. The 
correspondent free energy per site, — /3/q = jp. logQ | l 2 ^oo , follows in the form: 

1 r f dp dq 



~ l3fQ = 2J J 2^27 ln [( 1 + t ?)( 1 + t 2)-2ti(l-^)cosp-2t 2 (l-t 2 )cosgj, (9) 







while the true free energy per site, for Z = E exp(— f3H), can be recalculated from 
Z = (2cosh&i cosh6 2 ) i2 (5, and one finds: * 
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1 f f dp dq r I 

— f3f z = ln2-\ — / / In cosh 1b x cosh 26 2 — sinh2fe! cos p — sinh26 2 cos q , *(10) 

2 J J 2tt 2% l ' J 
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which is the famous Onsager's result, see Eq. (108) in [1]. An interesting comment by 
L. Onsager on the history of his remarkable solution can be found in [7]. The method we 
have applied, however, significantly differs from the original approach [1]. As it follows 
from the exact solution, in the ferromagnetic case, the critical point is given by the 
condition: 1 — 1\ — 1 2 — 1\ t 2 = 0, or sinh 2b± ■ sinh 2b 2 = 1. At this point, the (p — 0, q — 0) 
mode in (8)-(10) vanishes. 4 For the isotropic rectangular lattice, t\ = t 2 = t, the critical 

4 In the continuum-limit formulation, the parameter m = 1 — t\ — — t\ti plays the role of the 
Majorana mass of fermions [16,17]. It is interesting, that the Majorana-Dirac structures in 2D Ising 
model can be explicitly seen already at the lattice level, after the reduction of the integral (6) to two 
fermionic variables per site [16,17]. The zero mass corresponds to the critical point. 
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point is t c = tanh& c = \/2 — 1, or sinh2fe c = 1, with the inverse critical temperature 
b c = (J/kT c ) = \ ln(l + y/2) = 0.440686. The singularity in the specific heat near T c 
appears to be logarithmic, C/k ~ A c | logr| — > oo, with r oc |T — T c \ , r — > 0. For the 
isotropic lattice, the specific-heat critical amplitude is A c = ^ = 0.494539. For more 
comments on the properties of the 2DIM in a zero magnetic field see [1,6,17,18]. 

Let us now turn back to the factorized density matrix (5). The related identities, con- 
sidered below, may be of interest with respect to the 2D Ising model in a nonzero magnetic 
field, for which the analytic solution is yet unknown. Let L\ and L 2 be arbitrary linear 
forms in Grassmann variables, then: (1 + Li) (1+L 2 ) = exp (LiL 2 ) (1 + L 1 + L 2 ), where the 
nilpotent properties of fermions where taken into account. The two Grassmann factors 
here are combined into a one Grassmann factor accompanied by a Gaussian exponential, 
exp (LiL 2 ) = 1 + LiL 2 . The resulting identity can be iterated further on, and one finds: 

N —> N N N 

II (l + L a )= (l + £L Q ) exp (]T ]T L a L p ), * (11) 

a=l a=l a=l /3=a+l 

where L a are arbitrary linear forms in Grassmann variables, they also can include the 
spin variables as parameters [17]. In (11) we assume L N+ i = 0, similar conventions are 
assumed in (12)-(17) below. It is of interest also to consider the two mirror-ordered prod- 
ucts of factors like (11). They can be combined as follows: 



a a 
N > N < N 

n ( 1 + ^ )-n( i +^ 2) ) = ( i +E + ^) ) 

a=l a=l a=l 

N N N 



x exp { E E (4?> " ^) + Lf) + £ L^Lf } . (12) 

a=l /3=a+l a=l 

Making use of (11)-(12), the factorized density matrix (5) can be elaborated into a mixed 
spin-fermion Gaussian integral yet before the averaging over the spin variables. First, we 
apply the rule (11) to convert the local product of three factors A^B^A^ from (5) 
into a one Grassmann factor and accompanying Gaussian exponential. The exponensial 
factor in fact appears to be a part of the local exponential factor from (6) arising by 
the averaging over a mn at the junction in (5). The spin variables disappear from the 
exponential because of <r^ n = 1. Then we combine the arising two m-ordered products 
(at given n) according to (12) with respect to index m, and apply the rule (11) once again 
with respect to n. At the second step, the second line of (6) will be effectively completed. 
Finally, we come to the following representation for the spin-fermion density matrix (5): 

Q(a) = J D exp {S (a, b) + S mt (a, b,a)}, (13) 

where S (a,b) is the spin-independent part of the action, the same as in (6), so that 
the integral / Dexp (So (a, b)) is precisely the integral (6), while Si nt (a, b, a) is the spin- 
fermion part of the action, which explicit form is: 

L L L L L 

exp{S int ((7)} = exp{ E [ E L tnl L in<?mnCrpn + Yl E ^l^'^n^] } , (14) 

m=l n=l p=m+l p=l q=n+l 

whpre /~( 3 ) = (ft 1 ) n* -4- f( 2 ) h* A- n —b \ and T ( 4 ) = ffWn* _i_/( 2 )/i* _i_ 

" c c mn V mn m—ln ~ v mn u mn—\ ~ ^mn u mn! aLLK - L ^mn \ L mn m—ln ~ mn u mn—l ~ 

a mn + b mn )- The prefactor like 1 + L , where L is a linear form in fermions, is dropped 



4 



in (14), since, effectively, L = under the integral. It can be guessed that the averaging 
over the spins in (14) gives unity, as it is to be expected, assuming that there are no any 
other additional spin-dependent factors in the density matrix. The additional factors will 
appear in the case of a nonzero magnetic field. The expression (14) can be elaborated in 
such a way, that the spin variables can be eliminated as well in the presence of a field. This 
results, however, a theory with a nontrivial four-fermion integration, see (17). Schemati- 
cally, the spin-fermion part of action in (14) is of the form S = aaAa*a, where A is some 
nonlocal matrix. Introducing auxiliary fermionic fields at each site, c, c*, the action can 
be transformed as follows: S = aaAa*a — > cc* + a (aAc + c*a*) — > cA~ 1T c* + a (ac + c*a*). 
In the resulting action, the spin variables are coupled to fermions linearly, so the averaging 
over the spins is not a problem. Thus, introducing auxiliary fermionic fields, c mn , c^ n , 
after some transformations, for the spin-dependent part of action (13) we find: 

. L L L L 

exp {Si n t (c)} = / [| [| dc mn dc mn exp j ^ [ C mn C mn + ' mn C mn a mn + 

m=l n=l m=l n=l 
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+ a mn a mn ( C m+ln + • • • + C Ln ) ] + ^ C 2" ( C 0n+1 + • • • + Cq L ) j , CQ n = ^ C mn ■ (15) 

n=l m=l 

Introducing new variables of integration by c mn — > c' mn =c mn + c m+ \ n + . . . + c Ln , with 
the inverse transformation c mn = c mn —c m+ln , we obtain another form of the same identity: 

L L L L 

exp {S int (*)}=/ II II dc *nndc mn exp { ]T ]T [ c mn d rn c* mn + a mn (a rnn c m+ln + 

m=l n=l m=l n=l 
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+ C mn a mn) ] + C l« ( C l«+1 + • • • + C 1L ) } , d m C* nn = C* mn -C* n _ Vn . (16) 

n=l 
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In the case of a nonzero magnetic field the additional Boltzmann factors like + 
are to be introduced into the density matrix (5), where t^l = tanh/i m „, and h mn = (3H, 
is the dimensionless site-dependent magnetic field. This corresponds to adding the terms 
like (. . . + h mn <7 mn ) into the hamiltonian (1). The elimination of the spin variables in the 
presence of a field in (16) then results: 

. L L L L 

mn eX P { \_ C mn^m C mn ^ ^mn{ a mn C m+ln~^~ 

(°0 J m=l n=l m=l n=l 

L 

+ C mn a mn) + (^ — ^mn ) a mn a mn C m+ln C mn + C l" ( c lrt+l + • • • + C]x) |. (17) 

n=l 

Substituting (17) into (13), we obtain a purely fermionic expression for the partition 
function of the 2DIM in a nonzero magnetic field. The action that appears in (17) is non- 
Gaussian, we deal here with a system of interacting fermions. Therefore, the partition 
function can not be straightforwardly calculated. 5 There are few interesting features that 
can be observed in the representations like (14)- (17). We see that the long-range fermionic 
correlations emerge in the 2D Ising model in a nonzero magnetic field [17]. This statement 
can be formulated in few ways. In the identities like (14) or (15), we find explicitly the 
fermions of one sort to be coupled to the nonlocal sums of other fermionic variables (ac- 
companying spin variables may be assumed to be partly 'frozen' by switching on the field). 

5 For conformal field theory analysis of 2DIM in a magnetic field aiT~T c sec [14, f 5]. 
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In the representations like (16) and (17), this is elaborated into the form in which extra 
lattice fermions with zero mass (kinetic term c mn <9 m c,^ n ) are coupled to basic fermionic 
fields. Under the elimination of the spin variables in a nonzero field, on the other hand, a 
non-Gaussian quartic term arises, which prevents the exact solution. An approximation 
of the Hartree-Fock type has been applied to analyze the critical properties of 2DIM in a 
field in [17]. The singular part of the free energy near the critical isotherm, T = T c , h ^ 0, 
in the regime of 'strong' magnetic field, r 15 / 8 <K h «. 1, has been conjectured in the form: 



with m 2 oc r 2 , and A oc h M(r, h), where M(r, h) is magnetization, r oc \T C — T\ , r — > 0, 
and h is the magnetic field. The case A = exactly corresponds to the singular part of 
the free energy in a zero field that follows from the exact solution (9)-(10). Within such 
approximation, the singularity in the specific heat at the critical isotherm was found to 
be logarithmic, (C/k) sing — E c \ \nh | — > oo, h — > 0, with the amplitude E c = (8/15) A c , 
where A c is the thermal specific-heat amplitude at the critical isobar (T ^ T c , h = 0) 
[17]. For isotropic lattice, A c = (8/tt)6 2 = 0.494539, and E c = (64/15tt)6 2 = 0.263754, 
where b c = |ln(l + \/2). It might of interest to check these predictions by Monte- 
Carlo experiments. In physical aspect, it seems also to be highly desirable to clarify the 
mechanism of the spontaneous ordering in the 2D Ising model in terms of fermions. 

References 

1. L. Onsager, Phys. Rev. 65, 117 (1944). 

2. E. W. Montroll, R. B. Potts, and J. C. Ward, J. Math. Phys. 4, 308 (1963). 

3. T. D. Schultz, D. C. Mattis, and E. H. Lieb, Rev. Mod. Phys. 36, 856 (1964). 

4. L. P. Kadanoff and H. Ceva, Phys. Rev. B 3, 3918 (1971). 

5. T. T. Wu, B. M. McCoy, C. A. Tracy, and E. Barouch, Phys. Rev. B13, 316 (1976). 

6. K. Huang, Statistical Mechanics (Wiley, New York, 1987). 

7. L. Onsager, In: Critical Phenomena in Alloys, Magnets and Superconductors, 
R. E. Mills, E. Ascher, and R. I. Jaffe, editors (McGraw-Hill, New York, 1971), 
p. XIX-XXIV; p. 3-12. 

8. F. A. Berezin, Russ. Math. Surveys, 24, No. 3, 1 (1969). 

9. S. Samuel, J. Math. Phys. 21, 2806 (1980). 

10. C. Itzykson, Nucl. Phys. B 210 [FS6], 448 (1982). 

11. V. N. Plechko, Sov. Phys. Doklady, 30, 271 (1985). 

12. V. N. Plechko, Physica A, 152, 51 (1988). 

13. C. Itzykson and J.-M. Drouffe, Statistical Field Theory (Cambridge University 
Press, Cambridge, 1989). 

14. A. B. Zamolodchikov, Adv. Stud. Pure Math. 19, 641 (1989). 

15. G. Delfino and G. Mussardo, Nucl. Phys. B 455 (1995) 724. 

16. V. N. Plechko, Phys. Lett. A 239, 289 (1998). 

17. V. N. Plechko, J. Phys. Stud. (Ukr) 3, 312 (1999). A Special Issue Dedicated to 
90 th Anniversary of Professor N. N. Bogoliubov. 

18. B. M. McCoy, The 1999 Heineman Prize Address. Integrable Models in Statistical 
Mechanics: The Hidden Field With Unsolved Problems. - math-ph/9904003. 




(18) 



6 



